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ABSTRACT 

! 

During gas metal arc welding the metal within the heat affected zone 
undergoes a significant thermal transient which can result in changes of 
the microstructure of the material. Such changes can significantly affect 
the strength, ductility and resistance to corrosion of the material in 
this area. The object of this study was to develop and verify a model of 
the gas metal arc process which would predict the thermal transient such 
that the mechanical and metallurgical properties of the final base metal 
could be predicted. 

The equation of heat flow in two dimensions was modified by inclusion 
of terms to account for radiation and convection heat loss. In addition, 
the fact that material properties are a function of temperature was in- 
cluded. The resulting equation was highly nonlinear, but successfully 
approximated through use of finite difference approximations. The interior 
boundary of the temperature field was the boundary of the molten weld for 
each set of welding conditions. An empirical correlation was developed, 
therefore, to predict the size of the weld pool as a function of welding 
conditions. The model was verified by actually measuring the temperatures 
with thermocouples during welding. 

1 

j The predicted and measured temperatures were compared and found to be 
in excellent agreement whenever the weld bead was truly two dimensional. 
Peak temperatures were predicted within an accuracy of 5% or approximately 
132° F. 
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I. INTRODUCTION 



It was the objective of this thesis to develop a model of the gas 
metal arc welding process such that the temperature history of the metal 
could be predicted. The gas metal arc process was chosen because it is 
a type of arc welding that is sophisticated enough to be used to join 
complex, sensitive alloys. In addition, because high deposition rates 
are possible, and automatic welding machines are available, the process 
is compatible with production welding. 

Since simple models were already in existence which work well far 
from the weld centerline, the primary area of concern was the so-called 
heat affected zone. This is the part of the base plate which, although 
it does not become molten, has been heated to a high enough temperature 
to undergo significant detectable structural change. The final goal is 
to be able to predict the temperature history and therefore the cool-down 
rate and the final microstructure in the heat affected zone as a function 
of welding parameters. This information would be extremely useful in de- 
velopment of welding processes for new high strength steels such as 
HY-130. 

Development of the model was similar to that detailed in references 
5 and 7. First a pool size and shape correlation was developed by blowing 
out the weld puddles with high pressure argon and measuring the resulting 
holes. This information was plotted as a function of nondimensional 
parameters and an empirical relation derived. Second, the actual temper- 
atures were measured during welding operations. These values were compared 
with the results of the computer program. The input to this program 
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includes welding conditions and the size of the weld pool. The experi- 
mental and computed results were compared and found to be in reasonable 
agreement . 
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II. BACKGROUND 



The first real success in analytical prediction of temperature 
histories after welding was made by Rosenthal in the late 19 30 ’s. He ob- 
tained solutions for the governing heat flow equation for situations 
where the source is modeled as a point, line or plane in a quasi-steady 
state condition. Quasi-steady state here implies that the temperature 
distribution is measured relative to a moving coordinate system which 
has its origin directly under the electrode. His solutions were also 
subject to the following assumptions: 

(1) The physical and thermal properties of the medium are constant. 

(2) The radiative heat losses from the surface to the surroundings 
are small. 

(3) Any change of phase in the base metal and weld metal can be 
ignored. 

(4) As was mentioned, and most important of all, the heat source is 
modeled as a point, line or plane source. 

Experimental tests by Rosenthal and Schmerber in 1938 verified that 
the solutions found by Rosenthal were good sufficiently far from the weld 
centerline but broke down close to the weld since the molten puddle cer- 
tainly has dimensions. Rosenthal’s work is summarized in reference 6. 

In the period between 1940 and 1967 many investigators tried with 
varying degrees of success to improve on Rosenthal’s work. In 1952 
Wells (8) used the Rosenthal equations to predict the width of the weld 
bead. He used experimental work done by Jackson and Shrubstall to verify 
his conclusions. In 1955 Apps and Milner (2) did some work in confirmation 
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of the work done by Wells, but even more important for purposes of this 
work, they first expressed the theory that the total heat input to the 
base plate is composed of heat from the electric arc and heat from gase- 
ous conduction from the arc column. They were also the first to use high 
pressure gas to blow out the weld pools in an attempt to prove or disprove 
theoretical size and shape predictions. Wide variance between theory and 
experimental shapes was found. Adams (1) developed engineering relation- 
ships which reasonably accurately predict centerline cooling rate and peak 
temperatures in 1958. 

One of the major obstacles to the Rosenthal solutions is that all re- 
quire the heat input to the metal to be known. It is, of course, easy to 
measure voltage and current at the welding power supply, but, in addition, 
a value for the so-called welding efficiency must be estimated. In 1965 
Christensen et al. (3) in a detailed work measured arc efficiency. Un- 
fortunately the range of values for each welding process is very large. 

In addition in this work methods for finding length, width, depth, and 
cross sectional area of the weld bead were proposed. Christensen 1 s sug- 
gestions for nondimensional welding parameters are also important to this 
work. 

The importance of knowing the exact temperature histories in the heat 
affected zone was summarized nicely in the 1958 Adams lecture to the 
American Welding Society by E.F. Nippes. (4) Significant changes can oc- 
cur in the microstructure in this area. As metals become increasingly 
sophisticated and sensitive to high heat these changes become even more 
important. For very high strength materials welding processes must use 
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the optimal heat input to do the least amount of damage if the strength 
of the base plate is to be preserved. For these reasons the temperature 
history in the heat affected zone must be known in detail so that the 
cool-down rate can be computed and microstructure predicted through use 
of a continuous cooling transformation diagram. This idea was suggested 
by Nippes in 1958. 

The failings of all previous models have been that they are primarily 
based on the Rosenthal solution and its inherent approximations . As was 
mentioned previously these models are good far enough from the weld center- 
line, but become inaccurate in the heat affected zone which is the area of 
primary interest. It has been proven that to make accurate predictions 
the following must be considered: 

(1) The heat source in the metal is of finite size. 

(2) The effect of the variations of thermal properties of the 
material; 

(3) radiative and convective heat loss; 

(4) distributed heat. 

Distributed heat is that part of the energy which is transferred to the 
base plate through the gas cloud. The inclusion of these factors makes 
the original partial differential equation highly nonlinear. It has 
finally become possible to approximate the solution to such problems 
through the use of a large digital computer. 

The first truly accurate prediction of temperature in the heat af- 
fected zone was published in 1967 and 1968 by Tanbakuchi (7) and Pavelic. 

(5) Their method derives an empirical correlation for molten puddle size 
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and shape, and then uses the boundary of the molten pool as one boundary 
condition in a finite difference solution. Errors of less than 10% were 
achieved. The reason for their marked success was that not only were the 
factors listed above included, but their work was limited to the gas 
tungsten arc process, thus eliminating for the first time the dependence 
on an estimated welding machine efficiency. 

The principles developed by Pavelic and Tanbakuchi form the basis 
for this work. The gas metal arc process was chosen for previously men- 
tioned reasons. The spray transfer mode was chosen because it is com- 
patible with the higher welding speeds and with more different types of 
metals. Two shielding gases were tried, argon and argon with 5% oxygen. 
The latter was chosen to improve the stability of the arc and to reduce 
weld spatter. Low carbon steel was used because its properties are well 
known and it is readily available. The computer model was restricted to 
two dimensionality to reduce the computer memory required, which is very 
large even for the two-dimensional case. 
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III. POOL SIZE AND SHAPE CORRELATION 

Since any sort of finite difference algorithm requires accurate 
boundary conditions this aspect of the project was very important, and 
had as its goal the derivation of empirical relationships which would 
predict points on the melting isotherm of the weld puddle. Then, since 
the melting temperature of the weld and the ambient temperature are 
known, the required number of boundary conditions are known. A theoreti- 
cal analysis is stated in some detail in reference 7. The conclusion 
reached in his work was that an analytical method of determining the 
points of maximum length and width is possible. However, this analytic 
solution is tied to the 2D solution of the heat conduction equation and 
is not confirmed by experiment. As a result, in an attempt to profit from 
previous experience detailed in reference 5, the blow-out method of deter- 
mining size and shape of the puddle was tried first. 

A. Experimental Procedure 

All experimental work was done on a LINDE automatic gas metal arc 
system fed by a LINDE type SVI 300 welding power supply. After considera- 
ble trial and error a satisfactory blow-out system was constructed and is 
depicted schematically in Figure 1. The copper tubing was attached to the 
torch of the welding machine so that it had a constant orientation with 
respect to the arc and the molten puddle. The flow restricting valve was 
a coarse regulating needle valve and was left partially open. The blowing 
out of the puddle was accomplished by stopping the welding machine and 
opening the gate valve simultaneously. This apparatus was simple and un- 
sophisticated but was found to be effective once sufficient experience 
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Figure 1. PUDDLE BLOW-OUT SYSTEM 
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running the welding machine had been accumulated and once the correct 
angle of the final length of copper tubing was determined. The pressure 
of the argon used to blow the puddles out was not measured but is esti- 
mated to be 500 pounds per square inch. 

All the blow-out specimens measured two inches by 10.8 inches. This 
size specimen was chosen because it was narrow enough to be almost free 
of distortion during the welding and yet wide enough to prevent heat re- 
flected at the edges from affecting the welding conditions. Experience 
also showed the particular welding machine used consistently produced a 
stable weld after an initial rim of three to four inches. Therefore a 
length of over ten inches allowed some choice of the exact moment when 
the puddle was blown out. 

The following table summarizes the welding .conditions for the blow- 
out specimens. 



Voltage 


15-20 volts 


Current 


190-280 amps 


Speed 


28.5 to 48.7 inches per minute 


Thickness 


1/16 th inch 


Material 


mild steel 



The driving parameter is the speed. It was found that at each new -speed 
setting the voltage and current must be adjusted very exactly to get a 
two-dimensional puddle and yet to prevent the puddle from dripping through 
Attempts were also made to get two-dimensional puddles in l/8th inch steel 
but the increased width of the weld puddle allowed the molten material to 
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drip through before it could be blown out. This result is in agreement 
with Pavelic’s. Experimental results are included in Appendix C. 

B. Empirical Correlation 

By manipulation of the line source solution to the heat conduction 
equation a nondimensional welding parameter can be derived and is de- 
fined as follows: 

_ 3.415 El 
N5? = 27mtH m 

Details are included in Appendix D. This parameter is a function of the 
material properties and of the welding parameters of primary interest, 
namely voltage and current. It is not so obvious that the above parameter 
provides a method of excluding welding efficiency from further considera- 
tion for this procedure. If the dimensions of the puddle can now be ex- 
pressed in terms of N*, then for a particular machine, the isotherm of 
melting can be determined and the remainder of the solution is a standard 
finite element procedure consisting of averaging the temperatures of sur- 
rounding nodes until the solution converges. 

There is one flaw in this reasoning, however, and that is that the 
effect of the welding speed has not been included. A coordinate system 



is defined as follows: 
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where w is the moving coordinate in the direction of the welding and y 
is the coordinate used to express distance from the weld centerline. 
These coordinates can be nondimensionalized through dimensional analysis 
as follows : 

y* = XVy 
w* = XVw 

and the effect of welding speed has been included. 

In a manner similar to that used by Pavelic the blown-out weld pud- 
dles were measured using a Nikon machinist’s microscope. These results 
were fed into the service program "puddle fit" which is included in Ap- 
pendix B. As can be seen from Figures 2 and 3 the data falls roughly 
into an area which might be described by a relationship of the form 
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Figure 2. MAXIMUM PUDDLE LENGTH CORRELATION 



- 20 - 



m 




Figure 3. MAXIMUM PUDDLE WIDTH CORRELATION 



where C is the intercept of a straight line on a log log plot and m is 
its slope. Program M puddle fit" performed this operation using a least 
squares approximation. The results of this program are the following 
relations : 

for maximum length 

= 1.048 N" 1 * 359 



for maximum width 
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IV. REVIEW OF THE COMPUTER MODEL 



This section is a review of the computer model and a discussion of 
the minor changes made to it. The basic equations are unchanged. That 
is the model retains correction terms for diffused heat, convection, 
radiation, and variable material properties. Changes were of a minor 
sort and were intended to, first of all, make the program run in a speed 
range higher than that considered previously and second, to reduce the 
actual run time of the entire program. 

A. Basic Equations Used 

If it is assumed there is no heat generation in an isotropic solid 
and that the thermal properties are a function of t empe rat ure , the dif- 
ferential equation of heat may be written 

V • (kVT) • pc (4-1) 

Ds 

for a standard Cartesian coordinate system. If a quasi-stationary state 
is defined in which temperature is not a function of time and the coordi 
nate system is measured relative to the moving arc this equation is 
changed to the following: 

V • (kT)=-pcV^ (4-2) 

dW 

where w is the coordinate in the direction of motion and is described by 
the following: 

w = x - Vs 

If a normalizing factor is defined 

X = — -, m in 



V 



(4-3) 
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where S is a constant and will be discussed in detail later. Then all 
dimensions can be normalized by dividing by x q . Similarly temperatures 
are normalized by dividing by the difference between melting temperature 
and ambient temperature , that is , 



T* = 



T - To 



Tm -To 

After normalization, rearrangement, and introduction of the two- 
dimensional approximation equation (4-2) becomes: 

9 2 T* . u . -u - . - - • — = 0 



CL? 

N> 

51- 


h v x o 


3T* _ 1 


3k (T*) 


7 > 
3T* 


2 

i 


r 1 
3T* 


2 


3y*2 


0(T*) 


3w* k ( T" ) 


9T* 


3w : * 

J 


T 







(4-4) 

However, as was clearly shown by the differences in reference 5 and 
7, not all the energy is transferred into the base plate through the molten 
pool. Some energy is transferred directly from the ionized shielding gas 
to the base plate. A function which approximates this distributed heat is 
given below: 



q D (r) = q ( 0 ) e 



-Cr 2 



where r 2 = w 2 + y^ 



and which when normalized becomes 



qJ w *> y*) = 



q( 0 ) 



-Cx 2 (w'" 2 + y*2) 



D ' ' (T m - T 0 ) e 
Losses due to convection and radiation can also be included and the result 



is : 



9 2 T“ 


3 2 T“ 
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x o 2 

4 . SJ. - 
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1 — 1 
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1 3k (T*) 



k ( T* ) 3T" 



f . \ o 

3T* 2 



k 3 W" 



f . N 

3 T S! 



[3y* 



(4-5) 



In accordance with standard procedures the partial differentials are 
now replaced by finite differences which are now defined. 



3T* * ( a i ,i - a i~l> j ) T * . + a i~l» j 



TV . , . - 



3w*“ a i,j a i-l,j 1 ’'* a i,j ( a i , j +a i_i } j ) i+ljj 



a i,j 



a i-l,j (a i,j" a i-l, i" 1 ’] 



3T* b i,j b i-l,j , b i,j-l 

** = b i , j b i, j-1 l5j b i,j< b i,j +b i,j-l> 



b i>M + 



3 2 T* 
3 w*2 






1,1 



a i,j (a i,j + 






i+ l,j 



+ 



£ «p # 
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l-l, D i,] i-l,D 



3 2 T" 
3 y “ ^ 



b i,j b i,j-l 



L i ,1 



hi,j (b i , j + b i,j-l } Tl ’^ +1 



+ 



b. 

l 



, 1-1 




+ b i,j-i> 



ft 

T 

1,1-1 



where the computational module_ is defined in Figure 4. 
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Substitution of these expressions into equations and rearranging 
yields the following: 
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B . Convergence and Uniqueness 

In reference 7 Tanbakuchi used the following criteria for grid spac- 
ing parameters in order to insure convergence : 

9a • 
za min 

h < 

max a Q 

If h maxl which is a nondimensional grid spacing parameter in the program 
itself* is less than twice the minimum value of diffusivity divided by 
another value of diffusivity, then the solution will converge and will 
be unique. If h max is arbitrarily set to 1, then 

a = 2a . 
o mm 

and convergence and uniqueness will be satisfied if the length normalizing 
factor is defined 

^min 

x = (4-7) 

o y 

However, it was found that in the speed range greater than approxi- 
mately 30 inches per minute the length normalizing factor became very 
small which caused the puddle dimension to become very large. As a result, 
computer limits for division by very small numbers were exceeded in sub- 
routines TYPFN or COLUMN because k bessel function values returned from 

o 

subroutine BESSEL were very small. Similar problems had to be corrected 
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throughout the program as a result of the increased welding speed. The 
solution to this problem was to increase the length normalizing factor, 
x , by replacing the 2 in equation (4-7) by an arbitrary factor (see 
equation 4-3). This factor can be set at a value higher than 2 if the 
grid spacing parameter is set proportionally lower than 2. For most 
runs S equaled 4 and h max equaled .5. 

Unfortunately, however, such careful work does not insure convergence. 
As is pointed out by Pavelic, equation (4-4) is highly nonlinear and there- 
fore a finite difference approximation need never converge. The experi- 
mental results are available, however, for ready comparison so that a 
diverging solution should be immediately obvious. 

C. Changes in the Computer Program 

As was mentioned at the beginning of this section, the changes made 
to the program were minor, but worth mentioning. The first change was in 
the definition of the normalization factor and has been discussed in some 
detail. The second change was the establishment of two holding arrays in 
which successive iterations of the temperature surface are stored. After 
each complete iteration each element in each array is compared with the 
corresponding element in the other. If every element is within a toler- 
ance which can be specified in the input deck, execution stops and results 
are printed out. This improvement in the program reduced run time on an 
IBM 360 model 155 from several minutes to approximately one-and-one-half 
minutes. However, all the changes made so little difference to the struc- 
ture of the program that the flow charts in reference 7 still apply. 
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V. EXPERIMENTAL MEASUREMENT OF TEMPERATURE HISTORY 

Experimental temperatures were measured using the same welding ap- 
paratus described in section 3. Specimens were either l/16th inch or 
l/8th inch cold rolled SAE 1018 steel and measured 17.5 inches long by 
eight inches wide. Six thermocouples made of 24-gauge chromel and alumel 
wire were placed in each specimen. Filler metal was Airco type 675 which 
is recommended for low-carbon steel, and argon plus five percent oxygen 
continued to be used as the shielding gas. A Honeywell type 1508 record- 
er was used to measure and present the output of the thermocouples. A 
schematic of the measurement circuit is presented in Figure 5. 

A. Thermocouple Emplacement 

Welding conditions were first established using scrap steel. Then 
program puddle fit was used to estimate the width of the weld bead. Then 
thermocouples were placed at varying distances from the weld centerline, 
starting just outside the predicted width of the weld bead. Since temper- 
atures in the base plate can be predicted accurately by the point source 
theory as close as 1/4 of an inch from the weld centerline, all thermo- 
couples on all three temperature measuring runs were placed at distances 
less than .25 inches from the weld centerline. 

Since 24-gauge wire is approximately .020 inches in diameter, -thermo- 
couple holes were drilled with a diameter of approximately .040 inches. 
After the holes were drilled the wire was faced through the plate and the 
ends of the two wires were melted together. The small bead thus formed 
was hammered into the top surface of the plate to insure an intimate con- 
tact with the entire surface of the hole. After the thermocouples were in 
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Figure 5. 



THERMOCOUPLE CIRCUIT 
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place the actual distance from the weld centerline to the centerline was 
measured using the microscope mentioned previously. Figure 6 depicts a 
typical thermocouple installation. 

B . Temperature Measurement Procedure 

Figure 7 depicts the welding machine during a temperature measurement 
run. Voltage and current were measured by meters on the power supply. 

Speed was measured by using a stop watch to time travel over a known dis- 
tance. To start a run the Honeywell recorder was started at a speed of 
one inch per second. It is a feature of this recorder that timing marks 
are automatically made on the paper. Next the welding machine was started, 
and the start point of the welding machine was noted on the recorder paper. 
The welding machine was stopped several inches after the last thermocouple 
was passed, and the recorder after all temperatures had decayed sufficient- 
ly. In addition, for program ’’experimental data,” the distance from the 
start of the weld to the first thermocouple and the distances between 
thermocouples were measured. 

Three experimental runs were made with varying degrees of success. 
Welding parameters for each are listed below in Table 1. 



TABLE 1. 



Rim 


1 


2 


3 


Voltage (volts) 


15 


16.5 


22 


Amperage (amps) 


190 


235 


325 


Thickness (in.) 


1/16 


1/16 


1/8 



Speed (ipm) 



28.5 



43.5 



28.8 
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Figure 6. THERMOCOUPLE INSTALLATION 
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C. Analysis of Experimental Data 

The output of the Honeywell recorder is in the form of deflections 
continuously displayed. A sample curve is reproduced in Figure 8. The 
deflection was measured every half second and converted to a temperature 
by use of a calibration curve. These temperatures for an experimental 
run are included in Appendix C. This data plus the distances already 
mentioned form the input to program "experimental data" which puts the 
experimental data in the same coordinate system as the computer output 
for comparison. A listing of program experimental data is contained in 
Appendix B. 
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Figure 8. REPRESENTATION OF THERMOCOUPLE OUTPUTS ON THE HONEYWELL RECORDER 
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VI. COMPARISON OF COMPUTER AND EXPERIMENTAL RESULTS 

Welding parameters for the three experimental runs are listed in 
Table 1 in section V. Figures nine through 15 are a comparison of the 
temperatures measured in runs 1 and 3 with those computed using the weld- 
ing parameters and the thermocouple location. For information the thermo- 
couple locations are included in Appendix C. 

If error is defined: 

Tc* - Te* 

ERROR = (100%) 

then the error ranges in magnitude from as little as five percent at the 
maximum separation for thermocouple four to as great as 18 percent at 
the maximum separation for thermocouple one. This represents a differ- 
ence in temperature ranging between 132° R. and 476° R. at the points of 
maximum separation. V/ith this particular definition of error, the error 
observed is somewhat higher than that claimed by Pavelic in his work. It 
must be remembered, however, that his definition of error is based on a 
comparison of the peak temperatures regardless of when they occur rela- 
tive to each other. Using his definition the magnitude of the observed 
error is decreased considerably. 

As Pavelic pointed out the ideal comparison between experimental and 
computed data would occur if the curves were to overlap. However, in 
practice this is difficult to accomplish. Since experimental results 
generally lag the computed results , the contact between the wall of the 
hole and the thermocouple is immediately suspect. In fact, however, the 
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problem may be in the numerous factors and variables which are a part of 
welding and yet are extremely hard to quantify. 

As was mentioned before, welding conditions were determined first by 
choice of speed. The speed forced an optimal choice of voltage and cur- 
rent since the desired goal was a two-dimensional weld bead. However, 
as the speed or thickness increased, it was observed that the surface 
of the molten puddle flattened out and the resulting bead departed rad- 
ically from true two dimensionality. Experimental runs one and two were 
sectioned and macroscopically etched to determine the shape of the weld 
bead cross section. The results are represented in Figure 16. Since the 
thermocouples were a discontinuation in the specimen surface, all thermo- 
couples came in contact with the molten metal and consequently all gave 
relatively high readings. Figure 15 is thermocouple three of run three.. 
Using either definition, the error for this thermocouple is approximately 
20 percent and is unacceptably high, and this was the best thermocouple 
from this run. 

The problem in these two runs is that too much power is being trans- 
ferred to the metal. All welding handbooks consulted advised reducing the 
current and the voltage to reduce spatter and to reduce the tendency of 
the weld bead to spread out at the top. Of course, if the power were re- 
duced, the penetration is reduced and two dimensionality is lost anyway. 
One other possible variable is the type of welding gas. Argon plus 5% 
oxygen was used, however, which is recommended for gas metal arc welding 
on low-carbon steel. 

The conclusion reached was that the method first used by Tanbakuchi 
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and Pavelic is applicable to gas metal arc welding in a 
mode. However, a lesser degree of success was achieved 
ing speeds. 



spray transfer 
at very high weld- 



Thermocouple 
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TWO-DIMENSIONAL WELD BEAD (SLOW SPEED) 




NON-TWO-DIMENSIONAL WELD BEAD (HIGH SPEED) 



Figure 16. 



VII. SUMMARY AND RECOMMENDATIONS 



This method of prediction of temperature history in the heat af- 
fected zone is certainly useful because it is accurate. However, be- 
cause of the constraint of two dimensionality it is, in the opinion of 
the author, of only limited usefulness. For any sort of temperature pre- 
diction process to have any practical usefulness it must predict the 
temperatures under the actual welding conditions. Actual welding condi- 
tions in this case imply a somewhat lesser heat input per unit time than 
is necessary to achieve two dimensionality. The obvious implication is 
that this method should be extended into three dimensions. Only specula- 
tion is possible, however, whether the practical aspects of this method 
could be applied in three dimensions. It is doubtful, for example, that 
the blow-out system would work as well as it did. Since the blow-out 
equipment is obligated to be angled down toward the weld pool, it might 
be impossible to determine the depth of the weld pool to the high degree 
of accuracy required. 

The obvious answer to these difficulties is to make the angle between 
the plate and the blow-out nozzle as small as possible and increase the 
pressure of the blow-out gas , but only actual attempts would determine if 
the method is feasible. 

Unfortunately also, even the two-dimensional version used for this 
analysis required 400,000 bytes of memory to analyze an area one inch by 
three inches. It is obvious that this requirement grows rapidly as either 
more area is considered or as the third dimension is included. The need 
to have a large sophisticated computer greatly restricts the practical 
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value of this analysis since most production facilities like shipyards 
or aircraft manufacturing plants do not have such an installation. Of 
course if peripheral devices are available, they could be used to augment 
in core memory, but only if the program went through a major revision. 

With the advantage of hindsight, a far more reasonable approach seems 
to be to model the molten pool as a shape composed of many point sources. 
These point sources would be analytically varied in magnitude through the 
value of the molten metal to give the proper shape and heat input to the 
metal. This idea is similar to the analytical sources and sinks used to 
model shapes in the potential flow of fluids. The advantage here, if 
such an idea is possible, is that the temperature could be computed at any 
point in the material without an enormous amount of computer memory being 
required since the t emperature at every point is simply a combination of 
the temperatures resulting from the action of each small source. It is 
obvious that the above statement assumes linearity which would have to be 
checked experimentally. However, for practical applications, it is the 
opinion of the author that a way other than finite difference solution must 
be found in order to reduce the computer memory requirements. 
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A. Weld Puddle Correlation Data 
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37. 8 


.600 


.216 


9.38 


3.40 


4.83 


Fair 


6 


15 


210 


37.8 


.600 


.200 


9.45 


3.15 


4.83 


Fair 


7 


18 


250 


48.7 


.682 


.207 


13.84 


4.20 


6.90 


Good 


8 


18 


240 


42.5 


.591 


.206 


Not Computed 


- - - - 


Very Poor 


9 


16 


260 


42.5 


.600 


.205 


10.63 


3.64 


6.38 


Poor 


10 


15 


200 


32.4 


.584 


.229 


7.89 


3.10 


4.60 


Poor 


11 


15 


200 


32.4 


.576 


.211 


7.77 


2.84 


4.60 


Poor 


12 


15 


205 


32.4 


.575 


.201 


7.76 


2.72 


4.72 


Fair 


13 


15 


210 


32.4 


.575 


.184 


7.76 


2.48 


4.83 


Poor 


14 


15 


210 


32.4 


.632 


.198 


8.53 


2.67 


4.83 


Poor 


15 


15 


190 


28.5 


.598 


.185 


7.10 


2.20 


4.37 


Fair 


16 


15 


190 


28.5 


.583 


.192 


6.92 


2.28 


4.37 


Tood 


17 


17 


280 


48.7 


.728 


.209 


14.77 


4.24 


7.30 


Good 


18 


17 


2 80 


48.7 


.736 


.205 


14.93 


4.16 


7.30 


Fair 


19 


16.5 


2 35 


42.5 


.741 


.213 


13.12 


3.78 


5.95 


Good 


20 


16.5 


235 


42.5 


.784 


.190 


13.87 


3.37 


5.95 


Good 


21 


16 


230 


40.5 


.755 


.210 


12.61 


3.51 


5.64 


Good 


22 


16 


230 


40.5 


.728 


.183 


12.16 


3.06 


5.64 


Good 


23 


15 


210 


37.8 


.621 


.154 


9.78 


2.43 


4.83 


Good 


24 


15 


210 


37.8 


.650 


.191 


10.24 


3.01 


4.83 


Fair 


25 


15 


210 


37.8 


.606 


.172 


9.54 


2.71 


4. 83 


Fair 


26 


15 


205 


32.4 


.684 


.180 


9.23 


2.42 


4.72 


Good 


27 


15 


205 


32.4 


.627 


.196 


8.46 


2.65 


4.72 


Fair 


28 


15 


190 


28.5 


.684 


.198 


8.12 


2.35 


4.37 


Good 
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B. Experimental Temperatures 



Run #1 

Distance from 

Local Origin (in.) Temperature (°F.) at Thermocouple 





1 


2 


3 


4 


5 


6 


0.0 


77 


77 


77 


77 


77 


77 


0.5 


510 


330 


120 


90 


92 


90 


1.0 


875 


730 


240 


205 


290 


170 


1.5 


1130 


1005 


610 


490 


605 


355 


2.0 


1285 


1160 


810 


690 


725 


605 


2.5 


1370 


1285 


960 


790 


840 


710 


3.0 


1410 


1310 


1065 


890 


920 


805 


3.5 


1410 


1370 


1095 


970 


970 


890 


4.0 


1380 


1355 


1175 


1010 


1005 


955 


4.5 


1355 


1335 


1190 


1050 


1030 


990 


5.0 


1330 


1305 


1190 


1070 


1045 


1018 


5.5 


1290 


1290 


1185 


1085 


1050 


1030 


6.0 


1270 


1260 


1175 


10 85 


1050 


1040 


6.5 


1240 


1240 


116 5 


1070 


1040 


1050 


7.0 


1210 


1205 


1145 


1065 


1030 


1055 


7.5 


1190 


1180 


1135 


1050 


1025 


1050 


8.0 


1170 


1150 


1125 


1040 


1020 


1040 


8.5 


1150 


1130 


1105 


1025 


1010 


1040 


9.0 


1120 


1115 


1095 


1010 


995 


10 30 


9.5 


1090 


1090 


1080 


1000 


990 


1025 


Distance from 


Weld C/L. (in. ) 


.127 


.139 


.153 


.181 


.190 


.225 


XCCHRT(l) 


9.74 


12.6 8 


15.78 


19.15 


22.18 


24.95 



XCCHRT(l) is the distance in inches from the origin of the run to 
the local origin of each curve. 
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APPENDIX D. 



ANALYTICAL DETAILS 
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A. Origin of the Nondimensional Welding Parameter, N* 

The Rosenthal solution for the 2-D case assuming a thin plate and a 
line source is : 

At XVw - 



T - = — 



"t 2irk 

where Q = 3.415 nEI (Btu/hr) 



K 0 (XVr) 



( D-l) 



an 



d r = /w 2 + y 2 ( f t . ) 



H m , the heat content of the material, is defined: 

H m = < T m - V 

which ignores the heat required for the pearlite-austenite transformation 
and the latent heat of fusion which are small effects and very difficult 
to measure accurately. By multiplying and dividing equation (D-l) by 
pc (T m - T ) the following results: 



T - T = n 
o 



3.415EI 1 pc T m ~ T o 



2v 



k P=(T m - T 0 ) 



■XVw 



Ko (XVr) (D-2) 



which reduces to: 



T = n 



3.415 El 



-w 



2Trt H m a 



K 0 ( r* ) 



when the nondimensional quantities defined in sections III and IV are in- 
cluded. If a nondimensional welding parameter is defined: 



N* = 



3.415 El 



27ra t H 



m 



then the Rosenthal solution in the nondimensional form is the following: 



T* = n N* 



-W' 



K 0 ( r “ ) 



The same parameter v;as used for this study since it combines plate 
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thermal properties, thickness, welding voltage, and welding current in a 



convenient form. 
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